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ABSTRACT 

Geometric morphometric methods have made possible the study of shape and size inde¬ 
pendently. Here we have used uniform and non-uniform shape variables to study shell 
shape in two sibling species of Nassarius from NW Spain. The study revealed that the 
main shape difference between N. nitidus and N. reticulatus is that the former is more 
globular than the latter, a characteristic that is concordant with visual inspection. This 
method allows the use of multivariate and univariate statistical analysis on shell shape, as 
well as facilitating the biological interpretation of the shape variables obtained, and may 
thus be a fundamental tool for similar studies in taxonomy and evolution. 


RESUMEN 

Los metodos de geometria morfometrica han hecho posible el estudio de la forma y el 
tamano independientemente. Aqul hemos usado las variables de cambio de forma uni¬ 
forme y no uniforme para estudiar la forma de la concha en dos especies gemelas de 
Nassarius del NO de Espana. El estudio revelo que la principal diferencia de forma entre 
N. nitidus and N. reticulatus es que la primera es mas globular que la segunda, una 
caracteristica que esta en concordancia con la impresion visual. Este metodo permite el 
uso de analisis estadistico univariable y multivariable sobre la forma de la concha, asi 
como facilitar la interpretacion biologica de las variables asi obtenidas, y por lo tanto 
puede resultar una herramienta de gran utilidad para similares estudios de taxonomia y 
evolucion. 
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INTRODUCTION 

Definition and measuring of shape is 
not an easy task. Traditional morpho¬ 
metries used distances and ratios to 
characterize morphology of the shell in 
many species (Janson and Sundberg, 
1983; Johannesson, 1986; Mill and 
Grahame, 1995; Cruz, RolAn-Alvarez 
and Garcia, 2001). However, such an 
approach has some well-known deficien¬ 
cies (reviewed in Bookstein, 1991). First, 
selection of variables is very arbitrary 
and it can importantly affect the results 
obtained. Second, direct distances are 
typically highly correlated with the 
overall size of an organism, making it 
impossible in practice to study size and 
shape separately. Finally, even if ratios 
(between informative distances) were 
used as an attempt to correct the effect of 
size, such an approach has its own statis¬ 
tical problems (see Sokal and Rohlf, 
1996:17). Therefore,it was evident that 
classical methods were rather imperfect 
to single out size and shape components 
from general morphology. 

It was Sir D'Arcy Wentworth 
Thompson who first proposed depicting 
changes in shape by using reference de¬ 
formations of an organism superim¬ 
posed on a grid (see Bookstein, 1991). 
However, a general approach to solve 
this measuring quantitatively does not 
appear until shape could be described by 
the coordinates of a set of well-defined 
points or landmarks (Bookstein, 1991). 
Such landmarks link the geometry of the 
organism, the mathematics of deforma¬ 
tion and biological inference. Thus, a set 
of coordinates of homologous points 
from a group of individuals that have 
been adequately superimposed, i.e. 
translated, rotated and scaled to match 
one another by some reasonable criteri¬ 
on, can be used to compare and contrast 
their shapes. Then, differences in shape 
among individuals will be directly char¬ 
acterized by differences in these coordi¬ 
nates. The main aim of this paper is to 
introduce the taxonomist and marine 
ecologist to the use of geometric mor¬ 
phometric techniques, as they have sup¬ 
posed a revolution in morphometric 


studies (Adams, Rohlf and Slice, 2004). 
Oth&r alternatives to classical morpho¬ 
metric analysis s (based on index and ra¬ 
tios) are also known (see MacLeod, 
1999; Richtsmeier, Deleon and Lele, 
2002). Recently, the application of geo¬ 
metric morphometric methods to the 
study of shell morphology in Littorina 
saxatilis has opened the possibility of em¬ 
ploying it as standard methodology for 
morphological studies in molluscs (Car- 
vajal-RodrIguez, Conde-PadIn and 
RolAn-Alvarez, 2005). Here, we apply 
this technique to a well known pair of 
sibling species of Nassarius to character¬ 
ize their differences in shell shape. 

There has been considerable confu¬ 
sion in the taxonomy of the genus Nas¬ 
sarius from Atlantic shores (as an exam¬ 
ple, see Roding, 1798; Bucquoy, 
Dautzenberg and Dollfus, 1882; Lo- 
CARD, 1887; Coen, 1933; Nobre, 1940; 
Monterosato, 1912). Although Nassa¬ 
rius reticulatus (Linnaeus, 1758) and N. 
nitidus (Jeffreys, 1867) have been consid¬ 
ered by some of the former authors as 
sibling species (Jeffreys, 1867), the taxo¬ 
nomic consensus up to the early 1990s 
(Cernohorsky, 1984; Fretter and Gra¬ 
ham, 1984; Sabelli, Giannuzzi-Savelli 
and Bedulli, 1990; Poppe and Goto, 
1991) was to consider only N. reticulatus 
as a valid species. RolAn and Luque 
(1994) were the first to show that both 
were valid species, based on the compar¬ 
ison of museum types, as well as on new 
morphological and behavioural differ¬ 
ences detected in laboratory conditions 
studying two sympatric populations. 
They found that natural populations of 
N. reticulatus live preferentially along 
open shores and bays, while N. nitidus 
lives nearly exclusively in bays at NW of 
Spain. They also found both species liv¬ 
ing in sympatry in a few localities. San- 
juan, Perez-Losada and RolAn (1997) 
demonstrated the taxonomic status of 
these species by detecting diagnostic ge¬ 
netic (allozymic) differences in sympatry. 
A similar result has been observed with 
mtDNA variation in the same popula¬ 
tions (R. Barreiro, com. pers.). 

Thus, although there are diagnostic 
differences between these species in 
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Figure 1. Photo of N. reticulatus (A) and N. nitidus (B) in the established position for obtaining 
the shell measurements from the 13 landmarks (LMs) used (see text). 

Figura 1. Fotografia de N. reticulatus (A) y N. nitidus (B) en la posicion establecida para obtener las 
medidas de las 13 coordenadas (LMs) usadas (ver texto). 


individual shell characteristics (shell 
colour, profile, external lip, columellar 
callus, etc.; Rolan and Luque, 1994), 
shell shape differences between these 
species have not been formally demon¬ 
strated yet. In fact, the shell shape of 
both species is apparently rather similar 
(see Figure 1). Our interest in this study 
was twofold. First, to characterize shell 
shape in these two sibling species in 
sympatry, interpreting them in biological 
terms if possible. Second, to provide a 
representative example for taxomists 
and marine ecologists, that helps to 
apply this technique to similar cases. 
Our data demonstrate that this tech¬ 
nique can be of general utility as a diag¬ 
nostic tool to differentiate sibling species 
just in terms of shell shape (see also 
Johnston, Tabachnick and Bookstein, 
1991; Carvajal-Rodriguez et al., 2005). 


MATERIAL AND METHODS 

Adult specimens of both N. reticula¬ 
tus and N. nitidus were obtained from 
San Simon (Rfa de Vigo, Galicia, Spain) 


in 1992 and 1996. In this locality each 
species lives in different habitats, but 
usually buried in sand, sometimes 
under algae. All specimens within an 
area of a few square metres were cap¬ 
tured following the methods described 
by RolAn and Luque (1994), and taken 
to the laboratory. Then a subsample of 
20 adult snails (larger than 20 mm) for 
each species (40 in total) was chosen for 
morphometric studies. 

Shells were examined using a Leica 
MZ12 stereoscopic microscope, and 
colour images were captured and digi¬ 
tized using a Leica digital ICA video 
camera and QWin Lite version 2.2 soft¬ 
ware, always with the specimens in the 
same position (with the axis of the shell 
on the y-axis and the aperture on the 
same plane as the objective (Fig. 1) 
(Carvajal-Rodriguez et al., 2005). 
Shell variables were obtained using 12 
landmarks (LM) as shown in Figure 1. 
LM1 represents the apex of the shell; 
LM2 and LM3 are placed on the right 
border of the profile of the shell on the 
suture located at the beginning of both 
the penultimate and the last whorl. 
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Table I. Shape variables, percentage of variance explained and mean values and standard deviations 
within species. ANOVA check if there is a significant difference between species. ANCOVA check 
if the differences between species for the same variables are maintained in spite of size differences. 
Tab la I. Variables de la forma, porcentaje de la varianza explicado y valores medio y desviaciones estan- 
dar entre especies. El test ANOVA comprueba si hay diferencias significativas entre especies. El test 
ANCOVA comprueba si las diferencias entre especies para una variable se mantienen pese a las diferen¬ 
cias de talla. 



U1 

U2 

RW1 

RW2 

RW3 

RW4 

RWS 

RW6 

Variance Explained 

79% 

20% 

56% 

10% 

7% 

6% 

4% 

4% 

N.reticulatusN\m 

-0.015 

-0.002 

-0.038 

-0.001 

-0.001 

-0.001 

-0.002 

-0.001 

(std) 

(0.01) 

(0.01) 

(0.01) 

(0.01) 

(0.01) 

(0.01) 

(0.01) 

(0.01) 

N. nitidusN\m 

0.015 

0.002 

0.038 

0.001 

0.001 

-0.001 

,002 

0.001 

(std) 

(0.02) 

(0.10) 

(0.03) 

(0.02) 

(0.02) 

(0.01) 

(0.01) 

(0.01) 

ANOVA 

52.3*** 

1.9“ 

126.8*" 

■* 0.1 " s 

0.2 ns 

0.2 ns 

1.1* 

0.2 ns 

ANCOVA (size) 

25.7*** 

O.l" 5 

ns is non 

65.2**’ 
significant, ** J 

k 0.1 " s 
*p<0.001 

0.1 ns 

0.1 ns 

1.5“ 

0.1 ns 


respectively; LM13 and LM12 are the 
same points on the left side, respec¬ 
tively; LM4 is on the lower suture near 
the end of the last whorl; LM5 is on the 
upper part of the internal profile of the 
aperture; LM6 is the most external posi¬ 
tion on the external part of the outer lip; 
LM7 is on the most internal left side of 
the siphonal channel; LM8 and LM9 are 
placed on the most external point of the 
siphonal channel aperture on the right 
and left sides, respectively; LM10 is the 
confluence of the external siphonal 
channel with the last whorl; and LM11 
is the most external point on the last 
whorl (the periphery) at the left profile 
of the shell. Further details of the 
methods used are explained in Carva¬ 
jal-Rodriguez etal. (2005). 

For each specimen, centroid size 
(estimating shell size) and uniform 
(affine transformations, for example 
deformations which leave the sides of a 
square parallel) and non uniform (non- 
affine) components of shell shape where 
obtained. The estimation of shell shape 
components was accomplished by align¬ 
ing the raw coordinates of the speci¬ 
mens using the procrustes generalized 
orthogonal method (GLS; see Rohlf 
and Slice, 1990), with the scaling 


option of a= 0 (see further details in 
Carvajal-Rodriguez et al., 2005). We 
will not present here formal definitions 
or descriptions for these shell variables. 
Such information is clearly given in the 
book of Zelditch, Swiderski, Sheets 
and Fink (2004). Here we will focus on 
the biological interpretation of these 
variables. 

Centroid size is thus one of the best 
measures of size of any object indepen¬ 
dently of its shape (Bookstein, 1991). 
The two uniform variables (U1 and U2) 
describe variation affecting all land¬ 
marks simultaneously, after being cor¬ 
rected for scale, rotation, and superim¬ 
posed at the same position (Zelditch et 
al., 2004). An interpretation in geometri¬ 
cal terms of the differences at the 
uniform components is, however, not so 
simple (see Zelditch, Swiderski, 
Sheets and Fink, 2004). In our experi¬ 
ence with snails U1 has been typically 
associated with compression/dilation 
deformations (Carvajal-RodrIguez et 
al., 2005; this study). On the other side, 
the local component of variation gives 
rise to different variables that can 
explain orthogonal non-uniform compo¬ 
nents of the shape. When those non- 
uniform shape components are pre- 
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Figure 2. Representation of the two main local shape components (RW1 and RW2) in all speci¬ 
mens studied. The shape component RW1 is nearly diagnostic between the two species in these 
populations (see Table I). 

Figura 2. Representation de los dos componentes locales principals de la forma (RW1 y RW2) en todos 
los ejemplares estudiados. El aspecto del componente RW1 es casi diagnostico entre las dos especies en 
estas poblaciones (ver Tabla I). 


o Nassarius nitidus 
0 •Nassarius reticulatus 


sented relative to the overall shape 
information of the sample, as when 
using principal component analysis, the 
relative warps are produced (see 
Zelditch et al., 2004). The relative 
warps can be used as classical morpho¬ 
metric variables in statatistical univari¬ 
ate and multivariable analysis and tests 
(Zelditch et al., 2004). Furthermore, 
there is a possibility to represent the 
geometric meaning of the relative 
warps, using the thin-plate spline repre¬ 
sentation method (see Zelditch et al., 
2004). This gives the researcher the pos¬ 
sibility of interpreting an a priory 
abstract shape component in biological 
terms as we will show in this study. 

To get the landmark data in practice, 
we used digitalized photos of the shells 
as in Figure 1, and using the software 
TPSDIG (http://life.bio.sunysb.edu/ 
morph/morphmet/tpsdig2w32.exe) we 
obtained the coordinates for the different 
landmarks in the whole data set. Statisti¬ 
cal calculations to obtain centroid size 
and uniform and non-uniform shape 
components from landmark data were 
performed with the program TPSRELW, 
developed by Rohlf (1998) and avail¬ 
able at http://life.bio.sunysb.edu/ 
morph/morphmet/tpsrelww32.exe. 


Classical parametric tests were com¬ 
puted by the SPSS/PC package version 
12 . 0 . 1 . 


RESULTS 

Shell variation was decomposed in 
both size (centroid size: CS) and shape 
components. The two uniform shape 
components (U1 and U2) explained 79% 
and 20%, respectively of the uniform 
components of shape (after correcting 
for rotation and size effects). The vari¬ 
ables from the local shape component, 
the relative warps, explained a decreas¬ 
ing amount of the local variation (see 
the main six RW in Table I). Table I 
shows the univariate differences 
between species for those variables. 
Such analysis, however, assumes that 
size and shape components are indepen¬ 
dent (they do not show allometry). In 
our case, however, size still explained a 
significant amount of variation in shape 
(64%) using multivariate regression 
(Wilks' Lambda = 0.36; F = 7.04, P < 
0.001). However, the former differences 
between species are still detected cor¬ 
recting for size differences using an 
analysis of covariance (Table I). In 


85 







Iberus , 24 (2), 2006 


summary, we have detected differences 
in shape components irrespective of the 
differences in size between our two 
species. This is very convenient as dif¬ 
ferences in size can depend consider¬ 
ably on specimens available in the 
sample, especially when there is large 
within species variation in age or size in 
the field. The only limitation of the 
method is that we reduce our conclu¬ 
sions to the larger age/size classes in the 
species. 

Finally we quantified the degree of 
separation between our two groups 
using shape variables by a step-wise 
discriminant analysis. We observed sig¬ 
nificant differences between species for 
the discriminant function (Canonical 
correlation = 0.93; Wilks 7 Lambda = 0.13; 
_ 2 = 74.4, d.f. = 3; P < 0.001). A cross val¬ 
idation statistical analysis was also 
accomplished to verify the efficiency of 
the discrimination. All N. reticulatus and 
almost all (except one) N. nitidus speci¬ 
mens were properly classified a posteri¬ 
ori (more than 95% of success). Three 
shape variables contributed to the dis¬ 
criminant function (RW1, RW5 and 
RW7). In fact, the most efficient variable 
for discrimination was RW1 (see Figure 
2). In Figure 3 the thin plate spline rep¬ 
resentation allows us to interpret in geo¬ 
metric terms the positive and negative 
deviations values for RW1 between the 
two groups. The results graphically 
support the visual differences in shell 
morphology that can be observed 
between two representative specimens 
of these species. The advantage of RW1 
is that we can use a formal statistical 
variable which summarizes the main 
differences in shape (after correcting for 
size and rotation effects) between the 
two species. 


DISCUSSION 

We have shown that geometric mor¬ 
phometric methods can be used to find 
significant differences in shape between 
two sibling species of Nassarius. In this 
case the technique succeeded in discrim¬ 
inating in practice specimens of these 


two species just in terms of shell shape. 
This could be fundamental if one is 
interested in doing manipulative experi¬ 
ments directly in the field with a pair of 
sibling species, when diagnostic mor¬ 
phological traits are no available. In that 
case, the two species could be theoreti¬ 
cally determined by geometric morpho¬ 
metric methods and marked alive for 
future recapturing and monitoring. In 
addition, this method can be a valuable 
general tool to study and characterize 
shell and morphological differences in 
shape, with the possibility of interpret¬ 
ing shape in biological and simple 
terms. In the populations studied here, 
the main shape differences can be sum¬ 
marised in one shape variable RW1 
(Figure 3). Briefly, the figure shows that 
N. nitidus has a more globular shell 
shape than N. reticulatus (notice that this 
is a similar interpretation to the uniform 
differences observed for Ul; see Table I). 
This result is very intuitive as can be 
observed when visually comparing this 
two species (Figure 1), showing the 
good statistical properties of this 
approach. A similar result has been 
found when studying different ecotypes 
of Littorina saxatilis (see Carvajal- 
Rodriguez et al., 2005; Hollander, 
Lindegarth and Johannesson, 2005). 
Carvajal-Rodriguez et al. (2005), 
compared the efficiency of geometric- 
morphometric methods with the use of 
classical distances and ratios, and 
observed that the employment of land¬ 
marks considerably increases the chance 
of finding significant differences in 
shape between two groups. Further¬ 
more, even if classical principal compo¬ 
nent analysis made on many distance 
measurements (obtained from landmark 
data) can show similar statistical power 
to detect differences in morphology 
between two groups, that method can 
not safely discriminate between size and 
shape. Moreover, the significant princi¬ 
pal components detected are difficult to 
interpret in biological terms. In 
summary this technique has the main 
advantage of being an easy way to 
interpret in biological terms the differ¬ 
ences between our two species. 
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Figure 3. Thin-plate spline representation, from the TPSRELW software (Rohlf, 1998) for 
showing the most extreme positive and negative deformations of the landmarks with respect to the 
reference configuration for the main local shape component (RW1). Negative deformations are 
characteristic of N reticulatus and positive ones of N. nitidus (see Figure 2). 

Figura 3. Representation de “gradillas de deformation” del programa TPSRELW (Rohlf, 1998) 
mostrando la interpretation de las desviaciones mds positivas y negativas del primer componente local de 
forma (RW1) respecto a la configuration de referenda. Las deformaciones negativas son caracteristicas 
de N. reticulatus y las positivas de N. nitidus (ver Figura 2). 


Some caution is needed, however. As 
expected from geometrical and mathe¬ 
matical theory, all shape components 
should be analysed or interpreted 
together (Zelditch et al., 2004). This 
statement has some practical foundation 
as well, because as with principal com¬ 
ponent analysis, the percentage of local 
variation explained by RWs is negatively 
correlated with the number of such rela¬ 
tive warp. For example, in our data set 
the RW1 explained more than 50 % of 
the overall local variation in shape, and 
this percentage decreased to 10% for 
RW2 or even to 4% for RW6. A common 
sense argument would suppose that any 
shape difference between our species 
would be biologically more relevant and 
statistically safer if it were based on 
those relative warps that explain the 
largest amount of shape variability. 

The alternative position would be to 
avoid the use of any single relative warp 
just because of theoretical considerations 
(Zelditch et al., 2004). Such limitation, 
if followed strictly would suppose the 
loss of the main advantage of this new 
method and would limit its future utility. 


We believe that individual variables of 
shell shape can be extraordinarily useful 
to understand in biological terms how 
species are responding in evolutionary 
terms to environmental changes. The 
researcher should check in practice the 
real statistical importance of its geomet¬ 
ric variables as well as check experimen¬ 
tally the possible errors of the technique 
(see Carvajal-Rodriguez et al., 2005) 
before attempting to try a more detailed, 
and more interesting, biological interpre¬ 
tation of shape differences. 

We would like this work to con¬ 
tribute to the application of this new 
methodology to classical taxonomic and 
ecological studies in marine and terres¬ 
trial molluscs. 
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